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ABSTRACT 


This  thesis  is  a  design  of  an  alternative  system  which  may  provide  heating  to  the  Naval 
PostgraHnate  School  swimming  pool.  Particularly,  it  is  a  solar  heating  /  photovoltaic 
system  designed  for  a  better  efficiency  and  less  cost  of  installation  and  maintenance. 
Principles  of  heat  transfer,  control  and  fluid  dynamics  theory  are  used  for  the  determination 
of  this  heating  system  elements.  The  feasibility  of  its  installation  and  use  is  analyzed. 
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I.  INTRODUCTION 


A.  GENERAL 

Nowadays,  high  energy  consumption  based  on  fuel  has  been  associated  with 
economic  development  and  a  high  quality  of  life.  With  time,  energy  users'  demand 
increases,  while  the  fuel  reserves  are  dwindling.  The  conflict  between  demand  and  supply 
causes  an  escalation  of  the  fuel  cost,  which  should  mandate  the  research,  development  and 
use  of  alternative  resources  of  energy.  We  should  develop  alternative  means  of  energy  and 
fuel  conservation. 

B.  OBJECTIVES 

As  an  example  of  a  heating  usage  needing  optimization,  we  will  analyze  and  design  a 
new  heating  system  for  the  swimming  pool  of  the  Naval  Postgraduate  School. 

C.  SYSTEM  DESCRIPTION 

The  existing  system,  using  fuel  to  heat  the  water  and  electric  power  to  drive  the  pump 
for  circulation  through  the  heat  exchanger  and  the  filters,  is  old,  wastes  fuel  and  is 
inefficient.  Presently,  the  boiler  plant  is  used  to  keep  the  pool  at  the  desired  temperature. 
New  solar  system  designs  here  offer  several  alternative  energy  schemes  with  various 
tradeoffs.  An  optimum  design  is  suggested. 

The  objective  is  the  design  and  modeling  of  a  system  using  a  minimum  area  of 
collectors  and  electric  power  to  provide  a  maximum  thermal  energy  to  the  pool.  Taking 
into  account  its  geometry  and  masonry  materials,  as  well  as  accurate  weather  data  of  the 
Monterey  area  taken  from  meteorology  statistics  (Table  1),  and  applying  heat  transfer 
principles,  a  total  heat  loss  is  calculated.  That  heat  loss  is  balanced  by  the  collectors  and  the 
solar  irradiation  providing  thermal  energy  to  the  water  and  maintaining  its  temperature.  The 
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heat  loss  due  to  conduction,  convection,  radiation  and  evaporation  is  calculated  for  each 
month  in  Chapter  II. 

Power  vs  cost  tradeoffs  are  made  for  low-power  DC  pumps  with  photovoltaic  and  a 
separate  filter  pump  system  and  standard  AC. 


TAB 


LE  1. 


WEATHER  DATA  FOR  THE  MONTEREY  AREA. 


MONTH 

MONTHLY  TEMPERA-  (degrees  F) 

TURE 

WIND 

SPEED 

(mph) 

RELATIVE 

HUMIDITY 

9r 

Average 

maximum 

Average 

minimum 

Mean 

Average 

Average 

January 

60.0 

43.0 

51.5 

4  (1989) 

72 

February 

61.8 

44.5 

53.1 

4(1989) 

72 

March 

61.8 

44.6 

53.2 

4  (1989) 

73 

April 

63.4 

45.5 

54.4 

5 (1989) 

75 

EHI 

64.4 

47.6 

56.0 

6(1989) 

76 

June 

67.0 

50.1 

58.5 

6(1989) 

77 

68.0 

51.5 

59.7 

5 (1989) 

79 

HiH 

69.0 

52.6 

60.8 

5 (1989) 

79 

September 

72.0 

52.8 

62.4 

4(1989) 

78 

October 

70.3 

50.9 

60.6 

5 (1989) 

76 

November 

65.4 

47.0 

56.2 

4 (1989) 

74 

December 

61.0 

43.7 

52.3 

4  (1989) 

73 
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II.  HEAT  TRANSFER  ANALYSIS 


A .  HEAT  TRANSFER  DUE  TO  CONDUCTION 

It  is  considered  that  the  heat  is  transferred  by  conduction  under  one-dimensional, 
steady-state  conditions.  The  one-dimensional  refers  to  the  fact  that  one  coordinate  is 
needed  to  describe  the  spatial  variation  of  the  temperature,  and  the  heat  transfer  occurs  in 
that  direction  only,  as  can  be  seen  in  Figure  1.  The  system  (ground-walls  of  the  swimming 
pool  water)  is  characterized  by  steady  conditions,  since  there  is  no  energy  generated  per 
unit  volume  of  the  medium  of  the  wall,  assuming  that  the  temperature  is  independent  of 
time.  The  heat  diffusion  equation  is. 


y  tz 


+  q  =  pc 


n 

p 


(2.1) 


where, 

q=energy  rate,  which  is  generated  per  unit  volume  of  the  medium 
i3T 

pCp-^  =  time  rate  of  change  of  the  internal  (thermal)  energy  of  the  medium 
per  unit  volume, 
can  be  expressed  as; 


dx  V  dxj 


=  0 


(2.2) 


The  general  solution  of  the  eq.(2.2)  is  ;  T(x)  -  CjX  +  C2,  with  T(())  =  T^  ,  and  T(L)  = 
Tg  2-  Applying  these  conditions  to  the  equation  we  have  ^  C2  =  Tg  ,  and  at  x  =  L,  Tg  2= 
C^L  -t-  C2  =  CjL  +  Tg  j,  that  means  Cj  =  (Tg  2  -  Tg  ,  )  /  L. 
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Figure  1.  Heat  transfer  through  a  plane  wall. 


Substituting  into  the  general  solution,  the  temperature  distribution  is  then 


It  is  evident  from  this  result  that  the  temperature  varies  linearly  with  x  and  its  rate  of 
change  with  respect  to  the  x  coordinate  is: 


—  =  /'T  -T  'i  — 
dx  (  s,2  ^s.ijL 


(2.3) 


Then  the  conduction  heat  transfer  rate  is  obtained  by  using  the  Fourier's  law  equation. 


1  »  d'  1 '  j  V  k  ,  rr..  , 

qx=  -kA-^=  -  (Tg  2“^S,  l)  “  ir’(^S.  l“^S.2)  (2.4) 
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In  our  case  having  a  series  composite  wall  (Figure  2),  the  temperature  distribution  is 
the  sum  of  the  distributions  for  each  meuium  separately.  Therefore,  the  temperature 
distributions  for  the  insulating  fiber  resin  and  the  concrete  block  walls  are  T^fx)  and  T^,(x) 

respectively. 


T^(x)=:(Tj-Tp)^+Tp 


and 

T,(x)=(Tg-T,)j^  +  T^ 


(2.6) 


where,  Tj  is  the  temperature  at  the  thermal  contact  between  the  two  walls. 

Taking  the  thermal  conductivity  coefficients  of  the  fiber  resin  (k,-  0.038  W/mK)  and 
the  concrete  block  (kj,=  0.67  W/mK)  from  the  Tables  A-1  and  A-2  respectively,  the 
thicknesses  of  the  walls  (Lf=0.01m  and  L^=1.524m),  the  temperatures  of  the  water,  and 
the  ground  (T  and  T  )  and  applying  into  the  eq.  (2.4),  the  inlet  and  outlet  trtinsl'er  rates  are 

»  o 

obtained. 


'  X ,  in 


k^A 


and 


k,A 

lx,  out “  L  ^^1  " 


(2.7) 


(2.8) 
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Figure  2.  Temperature  distribution  into  the  swimming  pool  walls 


The  inlet  heat  rate  which  is  transferred  through  the  composite  wall,  from  the  hot 
water  (swimming  pool  water  at  80°F=300°K)  towards  the  co/d  ground,  which  is  assumed 
to  have  a  temperature  2°K  lower  than  the  mean  temperature  of  the  air,  is  equal  to  the  outlet 
heat  rate 


kfA 


kcA 


qx..„=qx.out «  — (Tp-T,)=— (T,  -Tg)« 


k,(T,-T,)  k,(T, -T^) 


(2.9) 
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Plugging  in  the  values  in  the  above  equation  we  obtain  the  temperature  Tj  and  then 
from  the  eq.  (2.7)  or  (2.8),  the  conduction  heat  transfer  rate  q^  is  calculated  for  each  wall 
and  the  bottom.  Summing  up  the  heat  transfer  rates  for  each  particular  area  A-,  the  total 
heat  loss  due  to  conduction  is  obtained  for  January  (Table  2).  Similarly  working  for  the 
mean  temperatures  in  Table  1,  the  conduction  heat  losses  for  each  month  are  calculated,  as 
shown  in  Table  3. 

_ TABLE  2.  CONDUCTION  HEAT  LOSS  (JANUARY). 


Aj  (22m^) 

(65m2) 

A,5  (652m2) 

A,^,,(849m2) 

qi(W) 

331.05 

478.19 

161.85 

478.19 

4796.64 

6245.93 

Therefore,  q^ond"  6245.93  Watt  =  15.86  MBtu  per  month  (for  January). 
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TABLE  3.  MONTHLY  CONDUCTION  HEAT  LOSSES. 


MONTH 

Watts 

MBtu  per  month 

January 

6245.93 

15.86 

February 

5780.91 

13.26 

March 

5764.18 

14.63 

April 

5540.04 

13.61 

May 

5242.26 

13.31 

June 

4777.28 

1 1.74 

July 

4553.14 

1  1.56 

August 

4349.06 

1  1.04 

September 

4051.32 

9.95 

October 

4385.86 

11.13 

November 

5205.5 

12.79 

December 

5931.46 

15.06 

B  .  HEAT  TRANSFER  DUE  TO  CONVECTION 

The  term  convection  refers  to  energy  transfer  which  occurs  between  a  fluid  in  motion  and  a 
bounding  surface  when  these  two  are  in  different  temperatures.  In  this  study  the  moving 
fluid  is  considered  to  be  the  atmospheric  wind  blowing  in  average  speeds  as  in  Table  1  for 
each  month,  and  the  bounding  surface  is  the  surface  of  the  water  of  the  pool.  Since  we 
have  a  flow  caused  by  external  means,  we  speak  of  forced  convection. 
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Considering  a  small  area  of  the  water  surface  with  temperature  T^,  the  wind  speed 

t* 

V,  and  air  temperature  X»;  a  local  heat  flux  q  occurs  which  is  proportional  to  the 
difference  (T^  -T„),  that  is: 

where,  h=local  convection  coefficient. 

•  i 

By  integrating  the  local  heat  flux  q  over  the  area  A^,  the  total  heat  transfer  rate  q  is 
obtained. 

q=lA 

^  S 

and 

q  =  (Ts-T.)JA  MA, 

_  • 

Then  defining  the  average  convection  coefficient  h  for  the  entire  surface,  the  total  heat 
transfer  rate  is  expressed  as: 

q  =  hA3(T,-TJ  (2.10) 

Setting  the  water  surface  area  A^=  A  =  652  m^,  the  water  temperature  T^,  =7^  =  80"F  = 
300°K  and  the  air  temperature  T  =  T^=  56®F  =  286.33°  K  (for  May)  the  eq.  (2.10) 
becomes: 

q  =  h  A(Tp-Ta)  (2.11) 
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— ►!  ^ -  dx  L 


Figure  3.  Heat  transfer  from  the  surface  of  the  water  to  the  air  due  to  the 

forced  convection. 


Assuming  that  the  wind  is  blowing  lengthwise  with  an  average  speed  as  shown  in 
Table  1,  we  calculate  the  Reynold's  number. 


R 


c 


VL 

V  . 
air 


where,  V  is  the  wind  speed,  L  the  length  of  the  swimming  pool,  and  the  kinematic 
viscosity  of  the  air.  From  the  Table  A-3,  for  the  mean  air  temperature  T^  of  each  month  we 
get  Vgjj.,  the  conductivity  coefficient  k  of  the  air  and  the  Prandtl  number  Pr.  Thereafter, 

calculating  the  Reynold's  numbers  for  each  case,  we  see  that  the  flow  over  the  water 
surface  is  turbulent,  because  5xl0^<Rg<10®.  Since  0.6<Pr<60  the  equation 


Nu=(0.037R^'’/5.g7i)pj.i/3 


(2.12) 
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is  used  [Ref.  3]  for  the  calculation  of  the  Nusselt  number  (Nu),  which  provides  a  measure 
of  the  convection  heat  transfer  occurring  at  the  surface.  The  Reynold  and  Prandtl  numbers 
are  applied  to  eq.  (2.12)  and  we  get  Nu.  As  described  in  [Ref.  3:pp.337],  the  defining 
expression  of  Nusselt  number  is: 


(2.13) 


Solving  eq.  (2.13)  with  respect  to  h’  and  substituting  the  obtained  values  of  k,  L,  and 
Nu,  we  get  the  convection  coefficient  h’.  Finally,  the  amount  of  heat  transfer  due  to 
convection  is  calculated  by  applying  the  values  of  h',  A,Tp,  and  to  the  eq.  (2. 1 1 ).  Then, 

we  have  the  results  for  each  month  as  shown  in  Table  4. 
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TABLE  4.  MONTHLY  CONVECTION  HEAT  LOSSES. 


MONTH 

Watts 

MBtu  per  month 

January 

50816.8 

129.01 

February 

47895.62 

109.82 

March 

47753.03 

121.23 

April 

55404.04 

136.11 

May 

60754.73 

154.23 

June 

54461.63 

133.79 

July 

43915.97 

1 1  1.49 

August 

41524.02 

105.41 

September 

31319.06 

76.94 

October 

41944.95 

106.48 

November 

42368.89 

104.09 

December 

49391.38 

125.39 

C.  HEAT  LOSS  DUE  TO  THERMAL  RADIATION 

The  heat  transfer  by  thermal  radiation  is  an  important  process  which  occurs  at  a  single 
surface,  in  this  situation,  and  depends  strongly  on  the  surface  geometry'  and  orientation,  as 
well  as  on  the  water  radiative  properties  and  temperature.  Thermal  radiation  is  related  to 
energy  released  as  a  result  of  oscillations  or  transitions  of  the  electrons  which  constitute  the 
water.  These  oscillations  are  sustained  by  the  internal  energy,  and  therefore,  the 
temperature  of  this  matter.  The  water  of  the  pool,  being  in  a  higher  temperature  T^  than 

that  of  the  surrounding  air,  tends  to  achieve  a  thermal  equilibrium  with  the  atmosphere  by 
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getting  cooler.  This  cooling  is  associated  with  a  reduction  of  the  internal  energy  stored  by 
it,  and  the  consequence  is  the  emission  of  thermal  radiation  from  the  surface.  This 
radiation  is  the  energy  rate  emitted  by  the  water  as  a  result  of  its  temperature  Tp  and  goes  on 

until  the  matter  reaches  the  atmosphere  temperature. 

The  nature  of  the  thermal  radiation  is  a  propagation  of  electromagnetic  waves  and  more 
particularly  those  with  a  wavelength  from  approximately  0.1  to  100  |im.  This  portion  of 
electromagnetic  radiation  spectrum  includes  a  part  of  ultra-violet,  all  of  the  visible  and  the 
infrared.  The  rate  at  which  the  radiation  leaves  the  surface  of  the  water  is  determined  by  the 
equation: 


q„d  =  AFae(Tp'- 


(2.14) 


as  in  [Ref.  3,  pp.  760],  where, 

A  =  Area  of  the  surface 
F  =  View  factor 
o  =  Stefan-Boltzmann  constant 
e  =  Emissivity  of  the  water  =  0.96  (from  Table  A-5) 

Tp  =  Water  temperature 

Atmosphere  temperature  =  273®K. 

The  view  factor  is  a  dimensionless  number  which  is  defined  as  the  fraction  of  the 
radiation  leaving  the  surface  of  the  water  and  is  intercepted  by  the  atmosphere.  For  this 
study  the  view  factor  is  considered  F  =  I.  However,  for  each  month  the  amount  of  heat 
transfer  by  thermal  radiation  is  the  same. 
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q„d  =  652m^  x  1  x  5.67  x  10"**  W  /  x  0.96  x  (SOO"*  -  273'*)  =  90336  Watt 

Then  the  monthly  heat  loss  due  to  radiation  is:  221.9  MBtu  per  month. 

D.  HEAT  TRANSFER  DUE  TO  EVAPORATION 

The  evaporation  from  a  free-water  surface  is  a  process  associated  with  water-vapor 
transfer  from  the  water  surface  to  the  air  and  occurs  by  the  combined  processes  of 
convection  and  diffusion.  Figure  4  is  the  schematic  model. 

The  surface  of  the  pool  at  temperature  Tp  is  exposed  to  a  stream  of  moist  air.  Into  the 
boundary  layer,  the  velocity  of  the  air  V(y)  varies  from  zero  to  at  the  outer  edge. 
Temperature  T(y)  decreases  from  Tp  to  T^,  and  the  humidity  decreases  in  the  y  direction 
from  (t>p  at  the  water  surface  to  in  the  air.  The  transfer  of  mass  of  water- vapor  having  an 
internal  energy  causes  a  thermal  loss  whose  rate  is  given  by  the  equation: 

^  cvap  ~  fg(  I  ■"  4*  a  )P  a,  sal  (2.15) 

where, 

hm  =  mass  transfer  coefficient  (m/s) 
hfg  =  heat  of  vaporization 
0^  =  relative  humidity 
Pa  jat  =  density  of  the  saturated  air 
A  =  area  of  the  swimming  pool. 


Figure  4.  Evaporation  of  water  into  moist  air. 


The  mass  transfer  coefficient  is  given  by  the  equation: 


h 


m 


_ h _ 

PaCp.a(Sc 


(2.16) 


where, 

h  =  convection  coefficient 
Pjj  =  air  density 
Cp  g  =  specific  heat  of  the  air 
Sc  =  Schmidt  number 
Pr  =  Prandtl  number. 


15 


The  values  of  the  above  variables  in  eq.  (2.15)  and  (2.16)  are  taken  from  [Ref.  3 
Tables  A-3  and  A-6]  and  correspond  to  the  mean  temperatures  of  each  month.  The  relative 
humidity  of  the  air  (J)^  is  the  mean  value  for  each  month  taken  from  Table  1 . 

However,  the  heat  transfer  due  to  evaporation  is  calculated  and  shown  in  Table  5. 


TABLE  5.  MONTHLY  EVAPORATION  HEAT  LOSSES. 


MONTH 

Watts 

MBtu  per  month 

18004.74 

45.71 

February 

19005.16 

43.58 

March 

18383.86 

46.67 

April 

21560.86 

52.97 

May 

25615.48 

65.03 

June 

27083.13 

66.54 

lljllllllflll 

22194.9 

56.34 

August 

23222.93 

58.95 

September 

21529.47 

52.89 

October 

26330.35 

66.84 

November 

19685.38 

48.36 

December 

17834.02 

45.27 
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III.  SYSTEM  OPERATION 


A.  SOLAR  INSOLATION 

The  average  solar  radiation  density  outside  the  atmosphere  is  approximately  1395 
W/m^  [Ref.  6]  for  a  plane  normal  to  the  rays  of  the  sun.  This  solar  constant,  reaching  the 
earth’s  surface,  reduces  because  of  the  atmospheric  reflectivity  and  absorption,  which 
depend  on  the  angle  of  the  incident  radiation  as  well  as  the  local  climatic  conditions.  In 
general,  the  average  reflectivity  is  35%  and  the  absorption  20%.  Then  the  average  incident 
power  density  is  approximately  S  =  (l-0.35)(  1-0.2)  1395  W/m^  =  725  W/m-  and  it  occurs 
at  noontime  (average  maximum  value).  Designing  a  solar  heating  system  is  of  interest  to 
calculate  the  daily  and  then  the  monthly  average  incident  power  density,  which  varies  with 
the  time  (seasonal  variation)  as  well  as  the  place  (geographic  position). 

1 .  Geometry  of  the  Incident  Solar  Radiation 

The  incident  solar  power  varies  due  to  the  daily  rotation  of  the  earth  and  to  the 
varying  tilt  of  its  axis  with  respect  to  the  equatorial  plane.  This  angle  varies  for  a  given 
place  in  a  sinusoidal  sense  (Figure  5)  through  the  year  and  gets  its  maximum  value  on  June 
22  (23.5®)  for  the  northern  hemisphere  and  its  minimum  one  on  December  22  (-23.5®). 
The  above  angle  is  called  declination  of  the  sun  0^.  Considering  a  spherical  coordinate 
system  (Figure  6)  with  its  origin  to  the  center  of  the  earth,  where  is  the  normal  vector  to 
the  surface  of  a  specific  place  M,  7t/2-6  is  its  latitude  and  is  its  hour  angle,  the 

incident  rays  of  a  flux  S  impinge  on  this  surface  with  an  angle  0^  with  respect  to  the 

equatorial  plane.  In  this  case,  for  convenience,  it  is  assumed  the  sun  rays  to  be  in  the  y  -  z 
plane,  so  noon  coincides  with  0  =  71/2. 
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Figure  5.  Declination  of  the  sun  through  the  year. 


Then  the  normal  incident  flux  is; 

S,  =  -a^S  (3.1) 

Resolving  the  vectors  a^.  and  S  into  rectangular  components,  we  have: 

sine  cos<l)  +  ay  sin  0sin<J)  +  a^cos0 
S  =  -ay  S  cosOj  -  a^  S  sin0^ 

and  substituting  in  eq.  (3.1),  we  get: 

Sj  =  S  cos0^j  sin0  sin(|)  +  S  sin0^  cos9  (3.2) 
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Figure  6.  Polar  coordinate  system  for  solar  radiation. 


The  above  equation  gives  the  incident  solar  radiation  with  respect  to  time  (0  and 
0j)  and  position  (tc/2-6).  To  obtain  the  average  daily  power  density  for  a  panicular  day  we 
integrate  it  over  the  daylight  time,  that  is: 


S 


r  av 


Ik 


j  (cos  6  j  sin  9  sin  4)  +  sin  cos  0)d0 

daylight 


(3.3) 


The  corresponding  hour  angles  for  the  sunrise  and  sunset  are  0^  and  (tt-0ni), 
respectively.  The  sunrise  angle  0^,  corresponds  to  that  0  which  makes  the  integrand  of  eq. 
(3.3)  zero,  that  is: 
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sin(J)m  =■ 

<}>m  =  sin 


-  1 


tan0 

r  tan  0 

Q 


tan  6  J 


(3.4) 


And  the  eq.  (3.3)  becomes; 


Srav  ~  **'”(cos0(isin0sin(j)  +  sin0(jcos0)d(j) 


S  r^/2,  „  ^  , 

=  — J  (cos0d  sin0sin0  +  sin0(j  cos0jd<J) 

TC  ^cn 

=  — [cos0d  sin0cos4)m  +sin0(i  cos0(;i  /  2  -  (j)^^)] 


(3.5) 


In  the  above  eq.  (3.4),  using  the  mentioned  solar  constant  (725  W/m^),  and 
substituting  0  =  53°  25'  (Monterey  latitude  =  36°  35'  N),  as  0^  the  mean  declination  of  the 

sun  for  each  month  taken  from  [Ref.  10]  and  as  daylight  period  the  mean  for  each  month, 
the  average  daily  solar  insolation  is  obtained,  as  shown  in  Table  6.  These  results  are 
consistent  to  the  data  based  on  weather  statistics  and  given  in  [Ref.  7]. 

2.  Optimization  of  the  Solar  Collector  Tilting 

In  designing  a  solar  heating  system,  the  solar  collector  tilting  is  taken  under 
consideration  in  order  to  maximize  the  received  solar  power  density.  The  basic  concept  for 
a  fixed  collector  is  that  it  must  receive  solar  energy  for  most  of  the  daylight. 
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TABLE  6.  MEAN  DAILY  SOLAR  INSOLATION  ( 

MONTEREY). 

Langley  /  day 

Btu  /  ft^  day 

January 

200 

737.4 

96.93 

February 

285 

1050.8 

138.12 

March 

430 

1585.41 

208.39 

April 

530 

1954.11 

256.85 

May 

575 

2120 

278.66 

June 

620 

2285.94 

300.47 

July 

600 

2212.2 

290.78 

540 

1991 

261.69 

September 

450 

1659.15 

218.08 

October 

350 

1290.45 

169.62 

November 

240 

884.88 

1 16.31 

December 

170 

627 

82.38 

Assuming  that  is  normal  to  the  collector  surface  and  a  is  the  angle  between 
the  normal  and  the  vertical  a^.,  then  the  incident  radiation  on  the  collector  is: 

S„=-a„S  (3^6) 

From  the  geometric  relations  (Figure  7)  substituting  in  eq.  (3.2)  the  angle  0  by 
(9  +  a),  we  get: 
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=  S  cos9j  sin(9+a)  sin0  +  S  sin9jjCos(9+a) 


(3.7) 


Integrating  eq.  (3.7)  over  the  daylight  period,  the  average  solar  irradiation 
is  obtained. 


^  [cos  9^  sin(9  +  a)sin  0  +  sin  9j  co.s(9  +  a)]d<|) 


daylight 


S 


~  TtT  i  [cos  9 j  sin  (9  +  a) sin  0  +  sin  9 ^  cos(9  +  a)]  d0 


=  — L  [cos  9 j  sin  (9  +  a)sin  0  +  sin  9 ^  cos(9  +  a)]d0 
=  -l^cos  9j  sin  (9  +  a)cos  0„  +  sin  9^  cos(9  +  a)  (^  -  0^)  (3.8) 
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In  order  to  find  the  optimal  tilting  a  of  the  collector  for  the  winter  (negative 
values  of  sun's  declination),  we  define  a  variable  which  is  the  sum  of  the  average 

solar  irradiation  for  each  day  over  the  winter  time  (Sep.  22  through  Mar.  22). 


Then 


S 


n  total 


X 

W  ini  er 


The  total  solar  irradiation  is  a  function  of  the  angle  a  and  a  specific  value  of  it 
maximizes  Numerically,  we  obtain  the  angle  a  after  a  series  of  runs  over  a  range 

of  feasible  angles  and  record  the  resulting  ,^,3,.  The  change  of  the  during  winter 
with  respect  to  the  tilting  angle  a,  is  shown  in  Figure  8.  It  can  be  seen  that  is 

maximized  for  a  =  56°. 


Figure  8.  Total  solar  irradiation  received  in  winter  (Sep.  22  through  Mar. 

22)  vs  tilting  angle  of  the  collector. 
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Tilting  angle  of  the  collector  (degrees) 


Figure  9.  Total  solar  irradiation  received  in  summer  (Mar.  22  through  Sep. 

22)  vs  tilting  angle  of  the  collector. 


If  the  collector  is  tilted  to  operate  during  the  summer  (Mar.  22  through  Sep.22) 
and  have  the  maximum  efficincy,  working  similarly  we  obtain  the  optimal  angle  a  =  10*^ 
(Figure  9).  For  a  year-round  operation,  the  optimal  angle  is  a  =  33®  (Figure  10). 
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Figure  10,  Total  solar  irradiation  received  by  the  collector  through  the  year 

vs  its  tilting  angle. 


B,  ENERGY  BALANCE 

In  Chapter  II  the  net  heat  rate  loss  due  to  heat  transfer  is  calculated  for  each  month 

(Table  7),  taking  under  consideration  the  weather  conditions,  the  geometry  and  the 
construction  of  the  swimming  pool.  Solar  irradiation  has  also  been  calculated  (Table  6) 

for  each  month  and  for  the  Monterey  area. 

Solar  power  is  collected  by  the  array  and  through  a  heat  exchanger  delivered  to  the 
thermal  load  (Figure  11).  Considering  the  average  efficiency  (fraction  of  the  solar 
irradiation  extracted  as  useful  energy)  of  a  collector  usually  used  in  this  type  of  heating 
systems  n  =  0.7,  then  the  power  drawn  by  the  collector  is  =  n  Q^.  Assuming  also  that 

the  heat  exchanger  provides  to  the  pool  all  the  thermal  energy  it  receives  from  the  array,  and 
additionally  its  free  surface  Ap  receives  directly  a  part  of  the  solar  power,  we  model  the 

collector  -  pool  system. 
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Figure  11.  Block  diagram  of  the  energy  balance  in  the  collector  -  pool 

system. 


Qc  +  ( 1  -pj)  Ap  -  Ql 


or 


n  Qg  +  (1-pg)  q^  Ap  — 


(3.9) 


where  Qg,  Q^,  are  energy  rates  (Watts),  Pj  is  the  water  reflectivity  taken  from  [Ref.  3], 
and  (l-pj)  is  the  fraction  of  the  solar  power  received  by  the  swimming  pool  surface. 

The  above  relation  denotes  the  energy  balance  that  must  exist  for  the  water  temperature 
to  be  maintained.  Eq.  (3.9)  becomes: 


n  Aj^q^  +  ( 1  "Ps)  qs  ~  Ql 


(3.10) 
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TABLE  7.  TOTAL  MONTHLY  HEAT  LOSS. 


MONTH 

NET  HEAT  LOSS 

(WATTS) 

MONTH 

NET  HEAT  LOSS 

(WATTS) 

January 

165403 

July 

161000 

February 

163018 

August 

159432 

March 

162237 

September 

147236 

April 

172841 

October 

162997 

May 

181948 

November 

157596 

June 

176658 

December 

163493 

where  q^,  Ql  are  power  densities  (W/m^),  Ps  =  3  and  A^^  is  the  area  of  the  collector  array 
(m^)  to  be  calculated.  Designing  the  heating  system  to  operate  all  year,  we  are  looking  for 
the  worst  case,  where  the  array  will  provide  energy  to  balance  the  greatest  loss  while  the 
incident  solar  irradiation  will  be  the  least  (maximum  apenure).  This  turns  to  be  -  2183 
m^  for  December.  If  the  heating  system  is  to  operate  from  April  through  October,  then 
=  720  m2  =  7750  ft^. 

The  above  calculated  areas  of  the  solar  collector  array,  which  is  to  be  used  in  each 
case,  are  too  large  and  off  the  swimming  pool  site.  Up  to  this  point,  the  design  of  the 
heating  system  was  done  without  taking  under  con.sideration  the  already  existing  swimming 
pool  cover.  Assuming  that  the  pool  will  be  open  from  10  a.m.  to  6  p.m.  and  covered  the 
rest  of  the  day,  then  the  net  heat  loss  Ql  will  be: 


Qi  =Q, 


cond.w 


1/3  Q, 


conv.srf 


2/3  Q, 


cond.cov 


1/3  Qevap  + 


1/3  Q 


rad.srf 


2/3  Q, 


ad.Cdv 


(3.11) 
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where,  Q^ond  w  conduction  heat  loss  through  the  walls  and  the  bottom  of  the  pool 
which  remains  the  same  after  covering  (Table  3),  Q^onv  srf  convection  heat  loss  from 
the  free  water  surface  which  holds  for  one  third  of  the  day  (open  pool)  (Table  4), 

is  the  conduction  heat  loss  through  the  pool  cover  which  holds  for  two  thirds  of  the  day 
(covered  pool),  is  the  evaporation  loss  for  open  pool  only  (Table  5),  is  the 

radiation  loss  through  the  free  water  surface  (one  third  of  the  day),  and  is  the 

radiation  loss  from  the  pool  cover.  This  cover  consists  of  two  polyethylene  films  and  a 
middle  cross  linked  polyethylene  foam  layer  (Figure  12). 

To  calculate  the  conduction  heat  loss  through  the  cover,  we  use  eq.  (2.7)  or  (2.8) 
taking  under  consideration  as  insulating  material  the  polyethylene  foam  middle  layer  only. 


Qconv.cov 


where  kp  is  the  thermal  conductivity  of  the  middle  layer  of  the  cover,  Lp  is  the  thickness  of 
the  layer,  Tp  the  temperature  of  the  pool  and  T^  the  average  temperature  of  the  air  and  A  the 

area  of  the  pool.  Taking  the  thennal  conductivity  coefficient  of  the  polyethylene  from  the 
Table  A-3  (k„  =  0.015  W/mK)  and  for  L  =  0.005  m,  T  =  300°  K,  T^  =  286.33°  K  (56° 
F  for  May)  and  A  =  652  m^,  we  obtain  Q^ondcov  “  46944  Watt  (for  May). 

For  the  calculation  of  Q^adcov’  (2- 14).  Then: 

Qrad.cov  =  AFae(Tp^-T^ini)  (3.12) 
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Figure  12.  Cross  section  of  the  swimming  pool  cover:  (a)  Polyethylene  scrim, 
(b)  cross-linked  polyethylene  foam,  (c)  polyethylene  film. 


where  A  is  the  area  of  the  swimming  pool  surface,  F  is  the  view  factor  and  F  =  1,  a  is  the 
Stefan-Boltzman  constant,  e  is  the  emissivity  of  the  outer  polyethylene  film  (£  =  0.85  taken 
from  Ref.  3),  =  300°  K,  and  273°  K  is  the  atmosphere  temperature.  Applying 

these  values  in  eq.  (3.12)  we  get  Qradcov  ”  79985.34  Watt. 

Substituting  in  eq.  (3.1 1)  all  the  values  of  the  different  losses  for  a  particular  month 
(e.g.,  for  May)  we  obtain: 


Ql  =  5242.26  -1-  0.33  x  60754.73  -i-  0.67  x  46944  +  0.33  x  2561 5.48  + 
+  0.33  X  90336  +  0.67  x  79985.34  =  148597.96  Watt. 


Comparing  the  above  result  with  the  total  heat  loss  for  May  from  Table  7  (181948 
Watt),  we  get  energy  savings  33350  Watt  or  18.32%.  Working  similarly  for  all  the  month 
heat  losses,  we  get  an  average  energy  savings  of  20%.  Then  using  eq.  (3.10),  we 
recalculate  the  area  of  the  collector  array  A^.  looking  again  for  the  worst  case  (greatest  loss  - 
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least  solar  irradiation),  while  tne  pool  cover  is  used  to  reduce  the  heat  loss  of  the  swimming 
pool  during  the  night  hours.  For  an  all-year  operation  of  the  heating  system,  it  is  required 
an  array  A^,  =  1616  m^.  If  the  system  operates  from  April  through  October,  then  A^  =  446 
m2  =  4800  ft2. 

For  a  high  performance  solar  heating  application,  as  in  this  case,  a  commonly-used 
collector  is  the  glazed  solar  collector  made  of  copper  plate  with  copper  piping,  (Figure  17). 
A  better  conductivity  of  the  thermal  energy  is  provided  by  copper  fins  roll-formed  around 
the  copper  tubes.  The  collector  surface  is  glazed  with  a  high-quality  glass  which  operates 
as  a  heat  trap,  reducing  radiative  and  convective  losses  emitted  by  the  collector.  The 
dimensions  of  a  formal  collector  panel  of  this  kind  are  4'xlO’;  for  the  heating  system 
operating  in  the  summer,  120  panels  are  required  and  their  cost  is  $  1300/panel  or  $156,0(X) 
for  the  whole  collector  array  including  its  installation. 

C.  PUMP  AND  PIPING 

The  solar  heating  system  includes  two  loops  of  fluid  circulation.  The  collector  loop 
consists  of  the  solar  array,  the  heat  exchanger  and  a  pump  Pj  for  circulating  the  water 

through  the  collector  and  back  to  the  exchanger  (Figure  13).  The  second  loop  circulates  the 
water  from  the  exchanger  to  the  filters  and  then  to  the  pool  using  a  pump  P^.  This  purrp 

already  exists  in  the  old  heating  system.  It  has  a  flow  rate  of  750  gal/min  circulating  the 
pool  water  roughly  three  times  a  day  and  is  driven  by  standard  ac.  The  other  pump  P,  is  to 

be  determined  with  respect  to  power  (h.p.),  that  is,  we  have  to  determine  a  flow  rate  (m^/s) 
for  achieving  a  good  conduction  of  heat  from  the  collector  to  the  water  as  it  passes  through. 
On  the  other  hand,  a  large  flow  rate  through  the  exchanger  is  a  poor  practice  resulting  in 
excessive  pumping  and  cost.  Then  Pj  must  be  adequate  for  delivering  sufficient  water 

flow  for  good  solar  collector  performance  during  the  summer  months  where  the  incident 
solar  irradiation  is  at  its  highest  (300.47  W/m^  in  June). 
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Assuming  that  all  the  heat  is  transferred  from  the  heat  collection  fluid  to  the  swimming 
pool  water,  then  the  thermal  energy  in  the  exchanger  is: 


Ee  =  (m^CcjAT,  -(mpCp)ATp 


(3.13) 


where  m^  is  the  water  mass  into  the  collector  loop,  c^,  is  the  specific  heat  of  the  water  in  an 
average  temperature  330”K  in  this  loop  and  c^  =  4184  J/kg  °K,  AT^  is  the  decreasing 
temperature  in  the  exchanger,  mp  is  the  water  mass  into  the  swimming  pool  loop.  Cp  is  the 
specific  heat  of  the  pool  water  in  an  average  temperature  300°K  in  this  loop  and  Cp  =  4180 
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J/kg  °K,  and  ATp  is  the  increasing  temperature  in  the  pool.  Dividing  both  parts  of 
eq.(3. 13)  by  time,  we  get  the  delivered  power  by  the  exchanger  to  the  load. 


Qe  =  (mcCc)ATc  -(mpCp)ATp 


(3.14) 


where  the  ^he  mass  flow  rates  (Kg/s)  provided  by  and  P2  respectively.  Eq. 

(3.13)  shows  that  the  rate  of  change  of  the  heat  in  the  exchanger  is  determined  by  the  rate 
of  heat  addition  from  the  collector  less  the  rate  of  heat  removal  by  the  swimming  pool 
(thermal  load).  Since  all  the  heat  is  transferred  from  the  collector  to  the  exchanger  we  have: 


nA^qs  =  Qe  =  (ihccjATc  -(mpCp)ATp 


(3.15) 


Determining  a  control  strategy  where  the  collector  pump  Pj  is  turned  on  when  AT^  = 
Tj,  •  Tp  in  =  o*"  ATp  =  2®K  below  the  desired  average  temperature  of  300°K,  and 

taking  under  consideration  that  rhp  =  750  gal/min  =  47.18  kg/s  and  q^  =  300.47  W/m^  (for 
June),  we  substitute  in  eq.(3.15)  and  solving  for  ihg  ,we  get  a  mass  flow  rate  29.17  Kg/s 
or  470  gal/min. 

Assuming  that  the  piping  of  the  heating  system  has  a  200'  length  of  5"  pipe  and  that 
the  fittings  (90s,  tees,  45s,  couplers)  have  an  equivalent  length  of  1(X)'  of  pipe  ,  then  for  a 
flow  rate  of  500  gal/min,  we  get  from  Table  A-7  a  head  loss  3x4.61'  =  13.83'. 
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Then: 


Total  head  loss  in  pipe  and  fittings  :  13.83’ 

Height  from  ground  to  roof  :  10.00' 

Approximate  friction  loss  in  plumbing  :  60.00' 

Total  head  loss  .  83.83' 

With  the  operating  -  system  characteristics: 

a.  Total  head  H  =  84  ft 


b.  Pump  discharge  U*-  470  gal/min 

we  have  to  determine  the  pump  Pj  size,  revolutions  and  break  horsepower  for  water 
circulated  at  an  average  temperature  330®K.  In  Figure  21  [Ref.  1 1]  using  the  measured  - 
performance  curves  for  a  particular  model  of  a  centrifugal  water  pump  with  an  impeller 
diameter  D  =  35  in  operating  in  710  r/min,  we  get  head  H  =  187.5  ft,  discharge  U  =  19000 
gal/min,  and  a  break  horsepower  P  =  1000  bhp. 

Given  the  above  data  for  a  model  pump  and  correlating  it  with  our  system 
characteristics  we  may  calculate  the  break  horsepower  of  P,  which  will  be  of  the  same 

geometric  family  (centrifugal)  and  operate  at  an  homologous  point  (same  performance). 
This  correlation  can  be  done  by  using  the  similarity  rules  [Ref.  1 1  [,  that  is,  head  goes  as 
n^D^,  discharge  as  nD^,  and  power  as  pn^D^,  where  n  is  the  number  of  pump  revolutions 
and  p  the  water  density.  The  similarity  rules  are: 


H 


V  "  / 


D 


(3.16) 


U 


n 

n 


vDy 


(3.17) 
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(3.18) 


3 

IdJ 

Applying  the  values  H,  D,  n,  taken  from  the  curves,  and  H*,  U*  for  our  system,  in  the 
eqations  (3.16),  (3.17)  we  get  an  impeller  diameter  for  D*  =  6.73  in  and  n*  =  2472 

r/min.  The  data  of  the  model  pump  taken  from  figure  21  corresponds  to  a  water 
temperature  of  60°F  or  288‘^K  with  density  p  =  1000  kg/m^.  In  the  collector  loop  (135°F 
or  330°K)  the  water  average  density  is  p*=  984  kg/m^.  Substituting  p*,  n*,  D*,  P,  p,  n, 
D  in  eq.  (3.18)  we  obtain  P*=  11  bhp  or  8.13  kW.  This  pump  costs  approximately 
$3,000. 

D.  CONTROL  UNIT 

The  controls  of  any  solar  heating  system  can  be  considered  its  brain.  With  the  aid  of 
sensors,  it  can  determine  when  heat  is  available  forcoUection  and  make  the  decision  when 
to  Stan  and  stop  the  heat  collection  process.  The  control  systems  used  in  such  solar  heating 
systems  are  called  differential  temperature  controllers.  Basically  they  compare  the 
temperature  of  the  solar  panels  to  that  of  the  load,  (Figure  14).  When  the  temperature  of 
the  panel  exceeds  that  of  the  load  by  4®C,  then  a  signal  is  sent  by  the  unit,  and  the  actuating 
devices  are  turned  on.  The  collector-to-load  differential  at  which  the  unit  turns  on  is  called 
turn  on  delta  T.  The  collector-to-load  at  which  the  unit  turns  off  is  called  turn  off  delta  T. 
The  difference  between  the  turn  on  and  turn  off  delta  T  is  called  temperature  hysteresis. 

Power  is  supplied  to  the  controller  from  120  Vac  household  power.  The  120  Vac  is 
converted  to  a  lower  value  dc  voltage  and  is  supplied  to  the  collector  temperature  detection 
circuitry  and  collector  sensor.  The  collector  temperature  detection  circuitry  receives  the 
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Figure  14.  Block  diagram  of  a  basic  differential  temperature  controller. 


collector  temperature  as  an  electrical  resistance  that  changes  with  temperature.  As  a  result 
of  this  resistance  and  dc  power  supplied  to  the  circuit,  a  voltage  is  developed  within  the 
detection  circuit  that  is  fed  to  the  comparing  circuitry.  The  dc  power  from  the  power 
supply  is  also  supplied  to  the  load  temperature  detection  circuitry,  which  operates 
identically  to  the  collector  temperature  detection  circuitry.  The  two  voltages  fed  to  the 
comparing  circuitry  are  compared  constantly.  When  the  comparing  circuitry  sees  that  the 
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Figure  15.  Differential  temperature  controller  output  wiring. 


collector  temperature  is  greater  than  that  of  the  load  by  a  set  amount,  it  send  a  signal  to  the 
output  circuitry,  i'.ie  circuitry  will  energize  a  relay  (Figure  15),  whose  constants  can  be 
sent  with  whatever  power  is  required  to  actuate  pumps,  valves,  etc.  When  the  output  of  the 
DTC  is  closure  of  a  set  of  relay  constants,  the  wiring  of  the  output  is  as  shown  in  Figure 
15. 

E.  PHOTOVOLTAIC  ARRAY  DESIGN 

The  photovoltaic  or  solar  cell  converts  incoming  light  into  direct  current.  Figure  16 
illustrates  a  construction  of  a  typical  photovoltaic  cell  which  is  very  similar  to  a  silicon 
diode.  When  the  front  and  back  contacts  are  connected  to  an  electrical  load,  current  will 
flow  as  a  result  of  the  induced  voltage  and  dc  electric  energy  will  be  generated.  The  solar 
cells  are  assembled  into  modules  which  are  the  smallest  units  that  the  consumer  will 
consider.  The  voltage  of  a  module  is  determined  by  the  number  of  solar  cells  connected  in 
series  and  usually  each  cell  contributes  0.5  Volts. 
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TOP 

SEMICONDUCTOR 


Figure  16.  Photovoltaic  (solar)  cell  construction. 


The  grid  -  connected  photovoltaic  system,  which  is  considered  as  the  electric  power 
source  for  the  pump  P,  and  the  control  unit,  consists  of  an  array  and  power  conditioner 

interconnected  with  ac  loads.  The  power  conditiong  unit  (PCU)  controls  the  operation  of 
the  array,  staning  it  up  in  the  morning,  shutting  it  off  at  night,  and  maintaining  the  proper 
voltage  for  maximum  power  extraction.  It  also  converts  the  variable  -  voltage  dc  electrical 
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output  of  the  array  to  constant  voltage  ac  electric  energy  for  the  load.  During  a  normal 
operation,  the  power  conditioner  controls  the  array  voltage  to  track  automatically  maximum 
power  whithin  certain  limits.  When  array  current  exceeds  the  capability  of  the  power 
conditioner,  such  as  may  occur  on  a  clear  day,  the  unit  will  increase  the  array  voltage  to 
limit  the  excess  current.  When  the  array  power  falls  below  a  minimum  level  the  PCD  will 
stop  automatically  converting  energy  and  shutdown. 

The  array  consists  of  several  series  strings  of  modules  which  are  then  paralleled  to 
achieve  the  desired  power  level.  The  series  strings  are  designed  to  provide  a  voltage  near 
the  center  of  the  power  condotioner's  input  voltage  operating  range.  Since  the  array 
configuration  depends  on  PCU  properties,  the  design  process  begins  by  choosing  the 
power  conditioner  and  then  designing  the  array.  To  achieve  a  desired  output  the  relation 
between  the  rating  of  the  power  conditioner  and  the  array  size  is  determined  by  the 
efficiency  of  the  power  conditioner  and  the  average  solar  insolation.  Since  the  usual  PCU 
efficiency  is  74%,  then  a  dc  input  of  about  1.35  kW  is  needed  for  each  IkW^^  and  is  given 

by  an  array  of  12  m^  (130  ft^)  approximately.  In  order  to  provide  ac  electric  power  8.5  kW 
to  the  pump  Pj  and  the  control  unit,  it  is  required  an  array  with  output  dc  power 

8.5kW/0.74  =11.5  kW.  For  a  formal  module  the  manifacturer's  specifications  under 
standard  operating  conditions  are  the  following: 

a.  Voltage  12  V 

b.  Current  2.5  A 

c.  Power  30  W 

Then  the  required  area  for  the  array  to  provide  the  desued  ac  power  is: 

2 

Sarrav  =  8.5kW  X  12 - =  102m^ 

array  . 
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Figure  17.  Interconnection  of  the  photovoltaic  system  with  the  load. 


Finally  determining  the  series  /  parallel  arrangement  of  modules  taking  into  account  the 
above  specifications, we  can  calculate  the  number  of  modules.  It  is  given  by: 


Array  Power  _  1 1500 
Module  Power  30 


=  383.3  =  384  module,^ 


while  the  number  of  modules  per  source  circuit  (series  connection)  is  given  by: 
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Array  Voltage  120 

- - - - —  = - =  10  modules 

Module  Voltage  12 

Then,  we  determine  the  number  of  modules  in  the  array  as  390  with  39  modules  in 
parallel  connection  producing  a  peak  current  97.5  A  and  10  in  series  providing  the  nominal 
voltage  of  120  Vdc  to  the  pump  and  the  control  unit  (Figure  17). 

Given  the  specifications  of  this  particular  module,  the  average  cost  of  a  panel  with  a 
power  of  110  Watts  is  $  500,  that  is,  the  cost  of  the  array  is  $  52,000  including  the 
installation. 
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IV.  CONCLUSIONS 


Economic  consideration  of  the  designed  solar  heating  system  is  very  important  to 
examine  and  analyze  its  economic  merit  and  feasibility.  The  alternative  must  be  compared 
with  the  already  existing  fuel  heating  system  on  an  equivalent  economic  basis. 

Using  records  and  information  from  the  boiler  plant,  we  calculate  the  average  amount 
of  steam  used  to  heat  the  pool.  From  boiler  data,  4400  fp  of  gas  produce  4925  lbs  of 
steam,  so  1  ft^  corresponds  to  1.12  lbs  of  steam  or  1084  Btu.  Also  each  therm  costs  $0.90 
and  the  average  amount  of  steam  is  20,000  lbs  per  day.  To  calculate  the  cost  per  day,  we 
use  the  equation: 


Cost  =  (lbs  of  steam  used)  x 


1  ft^  gas  1084  Btu 

- X - T - X 

L12  lbs  steam  ft  gas 


1  therm 
100000  Btu 


S  0.90 


X - 

therm 


which  gives  a  daily  cost  of  $175  or  $31,5{K)  for  a  season  of  180  days. 

For  the  designed  system  the  capital  cost  will  be  $156,000  +  $52,000  +  $3,000  = 
$211,000.  Assuming  a  service  life  of  10  years  with  an  interest  i  =  5  %  and  using  the 
straight  -  line  depreciation  plus  average  interest  method,  we  approximate  the  true  equivalent 
annual  cost.  These  two  variables  are  obtained  from  the  equations: 


Straight  -  line  depreciation 


P-L 

i 


(4.1) 


Average  interest  =  (P  -  L) 


n  + 1 


+  Li 


(4.2) 
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where, 

P  =  Capital  cost  (P  =  $21 1,000) 

L  =  Estimated  salvage  price  (In  our  case  we  assume  it  will  be  50  %  of  the  capital  cost) 
i  =  Interest  rate 
n  =  Estimated  life  time 

The  annual  cost  is  the  depreciation  plus  the  average  interest.  From  eq.  (4.1),  we  get  a 
depreciation  of  $21,100  and  an  average  interest  $8,176.  Then  the  annual  cost  is  $29,276. 

However,  with  our  new  designed  system  we  are  achieving  a  benefit  of  roughly  $2,200 
per  year.  In  this  evaluation,  the  maintenance  cost  and  the  personnel  payments,  which  are 
much  lower  for  a  solar  heating  system,  were  not  taken  into  consideration.  Its  efficiency  is 
also  dependent  on  weather  conditions,  but  not  on  fuel  price  variations. 
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APPENDIX  A  (HEAT  TRANSFER  TARLES) 
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TABLE  A-2  TIIERMOPII YSICAL  PROPERIES  OF  MASONRY 

MATERIALS. 


TYPICAL  PROPERTIES  AT  300  K 


DESCRIPTION/COMPOSITION 

DENSITY, 

P 

(kg/m*) 

THERMAL 
CONDUCTIVITY,  k 
(W/m  K.) 

SPECIFIC 
HEAT,  c, 

a/H  ■  K) 

Ruilding  Boards 

Asbestos-cement  board 

1,920 

058 

Gypsum  or  plaster  board 

800 

017 

— 

Plywood 

545 

012 

1,215 

Slicathing,  regular  density 

290 

0055 

1,300 

Acoustic  tile 

290 

0.058 

1,340 

llardboard,  siding 

640 

0094 

1,170 

llardboard,  high  density 

1,010 

015 

1,380 

Particle  board,  low  density 

590 

0.078 

1,300 

Particle  board,  high  density 

1,000 

0.170 

1,300 

Woods 

Hardwoods  (oak,  maple) 

720 

016 

1,255 

Soltwoods  (fir,  pine) 

510 

0.12 

1,380 

Masonry  Materials 

Cement  mortar 

1,860 

0.72 

780 

Brick,  common 

1,920 

0.72 

835 

Brick,  face 

2,083 

1.3 

— 

Clay  tile,  hollow 

1  cell  deep,  10  cm  thick 

0  52 

_ 

3  cells  deep,  30  cm  thick 

— 

0.69 

— 

Concrete  block,  3  oval  cores 
sand/gravcl,  20  cm  thick 

_ 

1.0 

_ 

cinder  aggregate,  20  cm  thick 

— 

067 

— 

Concrete  block,  rectangular  core 

2  cores,  20  cm  thick,  16  kg 

1.1 

_ 

same  with  filled  cores 

— 

060 

— 

Plastering  Materials 

Cement  plaster,  sand  aggregate 

1,860 

072 

_ 

Gypsum  plaster,  sand  aggregate 

1,680 

022 

1,085 

Gypsum  plaster,  vermiculite 

720 

0.25 

— 

aggregate 
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TARLE  A-3  TIIERMOPIIYSICAL  PROPERIES  OF  THE  AIR. 


T 

(K) 

P 

(kg/m’) 

(kJ/kg  ■  K) 

)i  •  10’ 

(N  •  s/m’) 

►  10‘ 
(m’/s) 

k  ■  10’ 

(W/m  •  K) 

a  10‘ 
(m’/s) 

Pr 

Air 

3,5562 

1.032 

71  1 

2.00 

9.34 

2  54 

0  786 

150 

2  3364 

ion 

1034 

4  426 

138 

5  84 

0758 

200 

1  7458 

1  007 

132.5 

7.590 

18  1 

10.3 

0737 

250 

1.3947 

1.006 

159,6 

11  44 

22  3 

15  9 

0720 

300 

1.1614 

1.007 

184  6 

15  89 

26.3 

22  5 

0707 

350 

0.9950 

1.009 

2082 

2092 

300 

29  9 

0700 

400 

0.8711 

1  014 

230.1 

2641 

33,8 

383 

0690 

450 

0.7740 

1  021 

2507 

32.39 

37.3 

47  2 

0.686 

500 

06964 

1  030 

270.1 

38.79 

40  7 

56.7 

0684 

550 

0.6329 

1.040 

288.4 

45  57 

43.9 

66  7 

0,683 

600 

0  5804 

1  051 

305.8 

52  69 

46  9 

769 

0685 

650 

0  5356 

1  063 

322.5 

60  21 

49  7 

87.3 

0690 

700 

04975 

1  075 

338.8 

68  10 

524 

980 

0695 

750 

04643 

1.087 

354.6 

76.37 

549 

109 

0702 

800 

0  4354 

1.099 

369.8 

84.93 

57.3 

120 

0.709 

850 

0  4097 

1.110 

384.3 

93.80 

59,6 

131 

0  716 

900 

0  3868 

1.121 

398.1 

1029 

620 

143 

0720 

950 

03666 

1  131 

411.3 

1122 

64.3 

155 

0723 

1000 

0  3482 

1  141 

424.4 

121  9 

66.7 

168 

0  726 

1100 

0.3166 

1.159 

449.0 

141.8 

71.5 

195 

0,728 

1200 

0,2902 

1,175 

473.0 

162.9 

76.3 

224 

0728 

1300 

02679 

1.189 

496.0 

185.1 

82 

238 

0719 

1400 

02488 

1,207 

530 

213 

91 

303 

0,703 

1500 

02322 

1,230 

557 

240 

100 

350 

0685 

1600 

0.2177 

1  248 

584 

268 

106 

390 

0.688 

1700 

0.2049 

1  267 

611 

298 

113 

435 

0.685 

1800 

0  1935 

1.286 

637 

329 

120 

482 

0683 

1900 

0.1833 

1  307 

663 

362 

128 

534 

0677 

2000 

0.1741 

1.337 

689 

396 

137 

589 

0672 

2100 

01658 

1.372 

715 

431 

147 

646 

0.667 

2200 

0  1582 

1  417 

740 

468 

160 

714 

0655 

2300 

0  1513 

1.478 

766 

506 

175 

783 

0647 

2400 

01448 

1.558 

792 

547 

196 

869 

0  630 

2500 

0.1389 

1  665 

818 

589 

222 

960 

0613 

3000 

0  1135 

2  726 

955 

841 

486 

1570 

0.536 

Ammonia  (NH,) 

300  06894 

2  158 

101  5 

147 

24.7 

16  6 

0887 

320 

0.6448 

2.170 

109 

169 

27.2 

194 

0870 

340 

0.6059 

2192 

1165 

192 

29  3 

221 

0872 

360 

0.5716 

2.221 

124 

21.7 

31  6 

24,9 

0872 

380 

0.5410 

2.254 

131 

242 

340 

27.9 

0  869 

400 

0.5136 

2.287 

138 

26.9 

370 

31.5 

0.853 

420 

0.4888 

2.322 

M5 

297 

404 

356 

0  833 

440 

04664 

2.357 

152.5 

32.7 

43  5 

396 

0826 

TABLE  A-4  THERMOPHYSICAL  PROPERIES  OF  THE  SA  ITJRATED 

WATER. 


TABLE  A-5  EMISSIVU  Y  OF  NONMETALLIC  MA  I  EUIALS. 


TEMPERATURE  EMJSSfVITV 
DESCRIPTION /COMPOSITION  (K)  r 


Silicon  carbide 

(») 

600 

0.87 

1000 

087 

1500 

0.85 

Skin 

(M 

300 

0.95 

Snow 

(/') 

273 

0.82-0  90 

Soil 

('0 

300 

0.93-096 

Rocks 

(/’} 

300 

0  88-095 

Teflon 

O') 

300 

085 

400 

0.87 

500 

0  92 

Vcgclation 

(M 

300 

G.92-096 

Water 

(M 

300 

0.96 

"Adapted  from  Reference  1. 

^Adaplcd  from  References  1.  9.  20,  and  21, 
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TABLE  A-6  THERMAL  AND  VAPOR  DIFFUSIVU  Y  FOR  DRV  AND 

SATURATED  MOIST  AIR. 


Tciiipcralure 

Degree  of 

a 

D 

r 

Saturalioii 

alD 

0 

0.770 

0.855 

50 

1 

0.769 

0.901 

0.854 

0 

0.799 

0.854 

60 

1 

0.797 

0.936 

0.852 

0 

0.828 

0.853 

70 

1 

0.826 

0.971 

0.850 

0 

0.858 

0.852 

80 

1 

0.854 

1.007 

0.848 

0 

0.888 

0.851 

90 

1 

0.883 

1.044 

0.846 

0 

0.919 

0.850 

100 

1 

0.911 

1.081 

0.843 

0 

0.949 

0.848 

no 

I 

0.938 

1.119 

0.838 

0 

0.981 

0,848 

120 

1 

0.963 

1.157 

0.832 

0 

1.012 

0.846 

130 

1 

0.985 

1.196 

0.823 

0 

1.044 

0.845 

140 

1 

1.003 

1.235 

0.812 
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I  ABLL:  A-7  head  loss  and  water  VELOCn  iES  INK)  IMI’ES 


Vclocily  mi-asuri'd  in  Loss  in  fis-t  of  waU-r  hcail  jirr  100  fl.  of  |>i|K', 


'A*  Va"  I'  'Vr  'Vi'  2-  z'/j'  r  3Vj'  4' 


AW.V.  I'r/ 

Lois 

Lott 

Wei 

Lott 

I'et 

£.«ii  1  >/  Lott 

1  Vel  Lott 

I'ei 

Lots 

I'ei 

Lott 

IV/ 

Lott 

5'./ 

Loti 

} 

2  74 

6  72 

1  4S 

1  51 

4 

5  48 

24  2 

2  V7 

>  45 

1  79 

1  54 

1  no 

59 

75 

177 

6 

K  25 

>1  2 

4  45 
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APPENDIX  B  (SWIMMING  POOL  PLANS). 


Figure  18.  Three-dimension  plan  of  the  Naval  Postgraduale  School 

swimming  pool. 
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Figure  19.  Heating  system  of  the  swimming  pool. 
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pump. 
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APPENDIX  C  (OPTIMIZATION). 


echo  on 
clc 

%  It  plots  the  total  solar  irradiation  (W/sq.m)  vs  tilting  angle  of  the 

%collector  for  the  winter  time. 

x=[50, 5 1,52,53,54,55,56,57,58,59,60,61,62]; 

>'=[41685,41760,41823,41873,41910,41934,41946,41944,41930,41904.41864.41812,4 

1747]; 

plot(x,y);grid; 

xlabel('Tilting  angle  of  the  collector  (degrees)'); 
ylabel('Total  solar  irradiation  (W/sq.m)'); 


echo  on 
clc 

%  It  plots  the  total  solar  irradiation  (W/sq.m)  vs  tilting  angle  of  the 
%  collector  for  the  summer  time. 
x=[4,5,6,7,8,9, 10, 11,12,13,14,15,16]; 

y=[44234,443 1 3,44378,4443 1 ,44469,44494,44506,44504,44488,44459.444 1 6,44360,4 
4290]; 

pIot(x,y);grid; 

xlabeK’Tilting  angle  of  the  collector  (degrees)'); 
ylabelCTotal  solar  irradiation  (W/sq.m)'); 
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echo  on 


clc 

%  It  plots  the  total  solar  irradiation  (W/sq.m)  vs  tilting  angle  of  the  %  ^collector  for  all 
over  the  year. 

x=[27,28,29,30,3 1 ,32,33,34,35,36,37,38,39]; 

y=[78680,78806,78907,78985,79038,79067,79072,79053,79010,78943, 78852.78737,7 
8598]; 

plot(x,y);grid; 

ylabel('Total  solar  irradiation  (W/sq.m)’); 
xlabel('Tilting  angle  of  the  collector  (degrees)’); 
echo  on 
clc 

a=0.0698 13 17007; 
theta=0.9322967; 

Snav=0; 

Sntotal=0; 

S=725; 
for  x=80:264; 

thetad=0.41015237422*sin(2*pi*(x-80)./365); 

phi=asin(-tan(thetad)./tan(theta)); 

Snav=(S/pi)*(cos(thetad)*sin(theta+a)’^cos(phi)+sin(thetad)*cos(theta+a)’^((pi/2)- 

phi)); 

Sntotal=Sntotal+Snav; 

end 

Sntotal 
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